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ABSTRACT 

Infrared Dark Clouds (IRDCs) are cold, dense regions of giant molecular clouds that are opaque at 
wavelengths ~ 10 /im or more and thus appear dark against the diffuse Galactic background emission. 
They are thought to be the progenitors of massive stars and star clusters. We use 8 /^m imaging 
data from Spitzer GLIMPSE to make extinction maps of 10 IRDCs, selected to be relatively nearby 
and massive. The extinction mapping technique requires construction of a model of the Galactic 
IR background intensity behind the cloud, which is achieved by correcting for foreground emission 
and then interpolating from the surrounding regions. The correction for foreground emission can be 
quite large, up to ~ 50% for clouds at ~ 5 kpc distance, thus restricting the utility of this technique 
to relatively nearby clouds. We investigate three methods for the interpolation, finding systematic 
differences at about the 10% level, which, for fiducial dust models, corresponds to a mass surface 
density S = 0.013g cm -2 , above which we conclude this extinction mapping technique attains validity. 
We examine the probability distribution function of S in IRDCs. From a qualitative comparison 
with numerical simulations of astrophysical turbulence, many clouds appear to have relatively narrow 
distributions suggesting relatively low (< 5) Mach numbers and/or dynamically strong magnetic fields. 
Given cloud kinematic distances, we derive cloud masses. Rathborne, Jackson & Simon identified cores 
within the clouds and measured their masses via mm dust emission. For 43 cores, we compare these 
mass estimates with those derived from our extinction mapping, finding good agreement: typically 
factors of < 2 difference for individual cores and an average systematic offset of < 10% for the adopted 
fiducial assumptions of each method. We find tentative evidence for a systematic variation of these 
mass ratios as a function of core density, which is consistent with models of ice mantle formation on 
dust grains and subsequent grain growth by coagulation, and/or with a temperature decrease in the 
densest cores. 

Subject headings: ISM: clouds, dust, extinction — stars: formation 



1. INTRODUCTION 

Large fractions of stars form in clusters from the dens- 
est clumps within giant molecular clouds (GMCs) (Lada 
& Lada 2003). These regions are also responsible for the 
birth of essentially all massive stars (de Wit et al. 2005). 
It is possible that our own solar system formed in such 
a region near a massive star (Hester et al. 2004), so the 
process of massive star and star cluster formation may 
be directly involved in our own origins. Understanding 
the formation of star clusters is also important as a foun- 
dation for understanding global galactic star formation 
rates (Kennicutt 1998) and thus the evolution of galaxies. 

In spite of this importance, there are many gaps in 
our knowledge of how massive stars and star clusters are 
formed. For massive star formation there is a basic de- 
bate about whether the process is simply a scaled-up 
version of low-mass star formation from gas cores (Shu, 
Adams, Lizano 1987), albeit requiring the high pressures 
found in the centers of star-forming clumps (McKee & 
Tan 2003), or whether a qualitatively different mecha- 
nism is involved, such as stellar collisions (Bonnell, Bate, 
Zinnecker 1998) or competitive Bondi-Hoyle accretion 
(Bonnell et al. 2001; Bonnell, Vine, & Bate 2004). A 
prediction of the scenario involving formation from cores 
is the presence of massive, gravitationally-bound starless 
cores as initial conditions. 

For star cluster formation there is a debate about 



whether the protocluster (or star-forming clump) is in 
quasi-virial equilibrium (Tan, Krumholz, & McKee 2006) 
or is undergoing rapid global collapse (Elmegreen 2000, 
2007; Hartmann & Burkert 2007). This corresponds to 
a debate about the timescale of star cluster formation: 
does it take many or just a few free-fall times? 

To help resolve these issues we require knowledge about 
the initial conditions of massive star and star cluster for- 
mation. The star-forming clumps that are the sites of 
these processes have typically been identified from the 
radiative emission and activity of their young stars: e.g. 
radio emission from ultracompact H II regions created 
by massive stars or maser emission from hot, dense gas 
near massive protostars. Unfortunately by the time this 
activity signposts the region, it has typically also erased 
much of the memory of the initial conditions from the 
system. This is especially true of protostellar outflows, 
which have a mechanical power large enough to signifi- 
cantly stir the gas of the star-forming clump and halt any 
large scale gravitational collapse (Nakamura & Li 2007). 

The initial conditions for massive star and star clus- 
ter formation are expected to be dense, cold cores and 
clumps of gas. In some formation models there may be 
of order hundreds of solar masses of gas compressed to 
within a few tenths of a parsec. For a spherical cloud the 
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For reference, £ = lg cm 2 corresponds to 4800M©pc 2 , 
A H = 4.3 x 10 23 cm" 2 and A v ~ 230 mag, for local dif- 
fuse ISM dust properties (e.g. Draine 2003). The exten- 
sion of the extinction law into the mid-infrared (MIR) is 
somewhat controversial and uncertain (Lutz et al. 1996; 
Draine 2003; Indebetouw et al. 2005; Roman-Zuniga et 
al. 2007), but nevertheless, such column densities are ex- 
pected to correspond to several magnitudes of extinction 
at 8 /jm. 

Such Infrared Dark Clouds (IRDCs) have been iden- 
tified in images of the Galactic plane from the Infrared 
Space Observatory (ISO) (Perault et al. 1996), the Mid- 
course Space Experiment (MSX) (Egan ct al. 1998). 
Simon et al. (2006a) identified about 10,000 potential 
IRDCs from intensity contrast features in the MSX sur- 
vey. Simon et al. (2006b) identified the 13 CO emission 
from about 300 of the darkest of these in the Galactic 
Ring Survey (GRS) (Jackson ct al. 2006), thus deriving 
their kinematic distances. Rathborne, Jackson, & Simon 
(2006) surveyed the 1.2 mm dust continuum emission 
from 38 of these IRDCs, identifying 140 mm-emission 
cores within these clouds. 

The temperatures in IRDCs are measured to be < 20 K 
(Carey et al. 1998). High dcuteration fractions have been 
reported by Pillai et al. (2007). Under these conditions 
one expects high depletion of volatiles onto ice mantles 
of dust grains (Dalgarno & Lepp 1984). Larger molecu- 
lar line surveys of IRDCs have been carried out by, for 
example, Ragan et al. (2006) and Sakai et al. (2008). 

In this paper we present a method of extinction map- 
ping of IRDCs using the diffuse Galactic IR emission ob- 
served by Spitzer us the background source. This method 
complements that based on measuring the extinction to 
individual stars (e.g. Roman-Zuniga ct al. 2007), by 
providing a measurement of the extinction at the loca- 
tion of almost every pixel in the cloud image, including 
at very high column densities, thus allowing a detailed 
investigation of the cloud structure. 

2. IRDC SAMPLE SELECTION 

Considering Spitzer Galactic Legacy Mid-Plane Sur- 
vey Extraordinaire (GLIMPSE) (Benjamin et al. 2003) 
8 /im (i.e. Infrared Array Camera [IRAC] band 4) im- 
ages, we chose 10 IRDCs from the sample of Rathborne 
et al. (2006), selecting those that were relatively nearby, 
massive, dark (i.e. showing high contrast compared to 
the surrounding diffuse emission), and/or with relatively 
simple surrounding diffuse emission. The properties of 
these clouds are listed in Table 1. Apart from being 
relatively nearby, this sub-sample is in fact fairly repre- 
sentative of the full 38 cloud sample of Rathborne et al. 
(2006). 

Simon et al. (2006a) fit ellipses to each cloud based on 
MSX images. While these ellipses are often not particu- 
larly accurate descriptions of the IRDC shapes, we will 
utilize them as convenient measures of the approximate 
sizes and shapes of the clouds, especially for the small 
scale median filter method of estimating the background 
radiation (§3.2.2). 



The Spitzer telescope has an angular resolution (PSF 
FWHM) of about 2" at 8 /im, which corresponds to a 
linear scale of 0.029 pc for a cloud at a distance of 3 kpc. 
Note, GLIMPSE images are processed to a pixel scale 
corresponding to an angular resolution of 1.2". 

3. IRDC EXTINCTION MAPPING METHODS 

The extinction mapping technique requires knowing 
the intensity of radiation directed towards the observer 
at a location behind the cloud of interest, I Vi o, and just 
in front of the cloud, I v ,\. Then for negligible emission 
in the cloud and a simplified ID geometry, 



= I v ,q cxp(-T„), 



(2) 



where the optical depth r„ = k„£, where k v is the total 
opacity at frequency v per unit gas mass and £ is the 
gas mass surface density. 

To evaluate n v appropriate for the various intensities 
measured by the Spitzer Space Telescope, namely the 
4 IRAC bands and the Multiband Imaging Photometer 
(MIPS) 24 (Um band, we adopt a spectrum of the diffuse 
Galactic IR emission from Li & Draine (2001), the filter 
response functions of the IRAC and MIPS bands, and 
perform an intensity and filter response weighted average 
of the opacity of various dust models (Figure 1 & Table 
2). Uncertainties in the dust models include the extent 
to which ice mantles have formed on the grains and the 
extent to which the grains have undergone coagulation. 

The IR background in the wavelength range probed 
by the IRAC bands receives its greatest contribution 
from the diffuse ISM (transiently heated small grains) 
at Band 4, i.e ~ 8 /um, compared to that from back- 
ground stars. Individual stellar sources become much 
more prominent in the GLIMPSE images at the shorter 
wavelengths. Thus in this paper we restrict our analysis 
to Band 4 images, leaving analysis at other wavelengths 
and the wavelength dependence of extinction for a future 
study. 

In the 8 /im band, we estimate a range of dust opac- 
ities per gas mass of 6.3 — 11 cm 2 g^ 1 , and adopt k v = 
7.5 cm 2 g" 1 , which is formally closest to the model of 
Ossenkopf & Henning (1994) with thin ice mantles that 
have undergone coagulation for 10 5 yr at a density of 
n.H = 10 6 cm~ 3 (or approximately equivalent to 10 6 yr 
at tih = 10 5 cm~ 3 , etc.). 

We estimate via interpolation from the regions sur- 
rounding the particular IRDC of interest. These inter- 
polation methods, which necessarily involve an averaging 
over small scale spatial variations in I v fi, are described 
below in §3.2. We evaluate from the observed inten- 
sities derived from the cloud images. First we consider 
the effects of foreground dust emission. 

3.1. Correction for Foreground Dust Emission 

Our determinations of both I v fi and I Vt \ are potentially 
affected by foreground emission from hot dust. We esti- 
mate the size of this effect given the (kinematic) distance 
to the cloud and a model for the Galactic distribution of 
hot dust emission, assuming it is the same as the distri- 
bution of the Galactic surface density of OB associations 
(McKee & Williams 1997): 
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TABLE 1 
Infrared Dark Cloud Sample 3 - 



Cloud Name 


I 


6 


Distance 




e 


P.A. 


/foro 


SsMF b 


-^LMF 






(°) 


(°) 


(kpc) 


(pc) 




(°) 




(g cm -2 ) 


(M ) 


(M e ) 


A (G018. 82-00.28) 


18.822 


-0.285 


4.8 


10.4 


0.961 


74 


0.209 


0.0355 


6,700 


7,600 


B (G019.27+00.07) 


19.271 


0.074 


2.4 


2.71 


0.977 


88 


0.075 


0.0387 


930 


830 


C (G028.37+00.07) 


28.373 


0.076 


5.0 


15.4 


0.632 


78 


0.266 


0.0527 


27,000 


42,000 


D (G028.53-00.25) 


28.531 


-0.251 


5.7 


16.9 


0.968 


60 


0.327 


0.0418 


17,400 


27,000 


E (G028.67+00.13) 


28.677 


0.132 


5.1 


11.5 


0.960 


103 


0.276 


0.0543 


15,100 


19,400 


F (G034.43+00.24) 


34.437 


0.245 


3.7 


3.50 


0.926 


79 


0.193 


0.0371 


1,770 


1,670 


G (G034.77-00.55) 


34.771 


-0.557 


2.9 


3.06 


0.953 


95 


0.140 


0.0420 


1,050 


1,140 


H (G035.39-00.33) 


35.395 


-0.336 


2.9 


9.69 


0.951 


59 


0.142 


0.0262 


4,000 


6,800 


I (G038.95-00.47) 


38.952 


-0.475 


2.7 


3.73 


0.917 


64 


0.141 


0.0616 


880 


1,490 


J (G053.11+00.05) 


53.116 


0.054 


1.8 


0.755 


0.583 


50 


0.121 


0.0699 


108 


80 



a Coordinate names, Galactic coordinates, kinematic distances, effective radii (of equal area circles), eccentricities and position angles of 
fitted ellipses are from Simon et al. (2006a). 

b This estimate of mean mass surface density, used to normalize the distributions in Fig. 11, is the areal average of those pixels for which 
values of Esmf > are derived. Estimates of a mean mass surface density based on Msmf and R e g are typically much smaller because 
of the regions inside the clouds ellipse with Esmf < 0. 



TABLE 2 

Spitzer Telescope Band and Background-Weighted Dust Opacities Per Gas Mass (cm 2 g _1 ) 



Dust Model a 


IRAC Band 1 
3.5^im b 


IRAC Band 2 
4.5^m 


IRAC Band 3 
5.9^tm 


IRAC Band 4 
7.8^m 


MIPS Band 1 
23.0^m 


WD01 R v = 3.1 


10.02 


6.26 


4.16 


6.25 


4.50 


WD01 R v = 3.1 flat IR bkg c 


9.78 


6.20 


4.25 


7.71 


4.24 


WD01 R v = 5.5 


15.56 


10.24 


6.55 


8.27 


5.54 


WD01 R v = 5.5 case B 


14.23 


11.49 


9.25 


10.96 


6.01 


OH94 thin mantle, yr 


17.52 (12.69) 


10.68 (8.69) 


7.78 (7.06) 


6.26 (6.15) 


4.43 


OH94 thin mantle, 10 5 yr, 10 6 cm- 3 


22.41 (16.23) 


13.31 (10.83) 


9.44 (8.56) 


7.48 (7.34) 


6.23 



a References: WD01 - Weingartner & Draine (2001); OH94 - Osscnkopf & Henning (1994), opacities have been scaled from values in 
parentheses to include contribution from scattering. 
b Mean wavelengths weighted by filter response and background spectrum. 
c No weighting made for spectrum of Galactic diffuse emission. 



where R is the galactocentric radius and Hr = 3.5 kpc 
is the radial scale length. For each IRDC, given its dis- 
tance and Galactic longitude, we calculate the ratio of 
the column of hot dust between the solar position (at 
R = 8 kpc) and the total column extending out to a 
galactocentric radius of 16 kpc. This "Foreground Inten- 
sity Ratio" , /foro, is listed in Table 1 for each IRDC. Then 
we derive an estimate of the true intensity of the radia- 
tion field just behind the cloud, J^o, from that measured 
via interpolation of the cloud images, I v $, obs, via 

Ivfi = (1 — ,/forc)^,0,obs- (4) 

We also estimate the true intensity of the radiation field 
just in front of the cloud, I Vj \ from that measured directly 
from the cloud images, I^^obs, via 

Iv.l — -^[Al.obs — ./forc-^iAO.obs- (5) 

Typical values of /f ore are about 15%, with values up 
to 33% for the most distant cloud, D, at 5.7 kpc. As 
an example of the size of the corrections to £ result- 
ing from the foreground subtraction, consider cloud C, 
for which /f ore = 0.266, Ii,,o,obs — 100 MJy sr" 1 and 
^!/,i,obs — 50 MJy sr -1 in the darkest part of the cloud. 
For these regions, one would derive £ ~ 0.092 g cm~ 2 
without the foreground correction and £ ~ 0.152 gem -2 
with this correction. 



Our estimate of /foro is uncertain due to small scale 
spatial variations in the hot dust emission in the Galaxy. 
Also, we are typically measuring I„,o,obs from regions 
relatively close to the IRDC of interest. Such regions 
are likely to overlap with the GMC hosting the IRDC, 
and thus have higher than average extinction of the in- 
tegrated Galactic background emission. This will cause 
us to tend to underestimate /foro and thus £. An upper 
limit to /fore is provided by the minimum value of i^,i i0 bs 
for each cloud. For example, for Cloud C this is about 
40 MJy sr -1 intensity units, so that for a background of 
100 MJy sr -1 , the maximum value of /foro — 0.4. Uncer- 
tainties in /fore are one of the major reasons for restrict- 
ing the extinction mapping analysis to relatively nearby 
clouds. 

3.2. Background Estimation 
3.2.1. Large-Scale Median Filter (LMF) 

A relatively simple way of estimating the diffuse IR 
background at a given location behind an IRDC is to 
take a median average of a region (i.e. filter) centered 
on the location of interest and that is large compared 
to the cloud. This method was applied by Simon et al. 
(2006a) to model the Galactic background from MSX 
images to then identify IRDCs as high contrast features. 
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Fig. 1. — Top: Adopted spectrum of diffuse fR background from the "MIRS" region in the Galactic plane at I a 44° modeled by Li 
& Draine (2001). Middle: Filter response of IRAC bands 1-4 and the MIPS 24 fim band. Bottom: Dust opacities (per gas mass) from 
Weingartner & Draine (2001) as updated by Draine (2003): Model with Ry = 3.1 (solid line), Ry = 5.5 (dotted line), Ry = 5.5 case 
B (dashed line). These opacities include scattering, which is about a 40% of total effect at 3 (im, but only a 2% effect at 8 ^m. The 
long-dashed line shows the thin ice mantle MRN model (only contribution from ft a b s per gas mass is shown) of Ossenkopf & Henning (1994). 
The dot-dashed line shows this model after 10 5 yr of coagulation at a density of njj = 10 6 cm -3 . The background and filter-weighted 
opacities are listed in Table 2. 



This method will only capture background fluctuations 
on scales larger than the cloud. Indeed if the filter be- 
comes too small, i.e. with a size comparable to the cloud, 
then the derived background will become influenced by 
the cloud itself and will be underestimated. 

For simplicity, we adopt a square-shaped filter for this 
method. After some experimentation and given the sizes 
of the IRDCs in our sample, we chose a filter size of 13', 
i.e. 650 Spitzer-GLIMPSE image pixels. We evaluate 
the background on a uniform grid with spacings of 24". 
At each location, we first consider the full distribution 
of pixel intensities, and evaluate the mode. Since the 
high-end tail of this distribution, due to extended bright 
sources and stars, was typically much larger than the 



low-end tail, due to the IRDC, we excluded pixels that 
had intensities twice the value of the mode, and then 
found the median value of the remaining distribution. 
This value was assigned to 7„ ; o, obs for all the pixels in 
the 24" x 24" square sharing the same center as the 13' 
filter. This Large Scale Median Filter (LMF) method is 
illustrated for cloud C in the left-hand panels of Figure 2. 
The distribution of pixel intensities is shown in Figure 3. 

Note that there will be some regions where the esti- 
mated value of /„,i is greater than that of I„q. This 
can arise because of foreground or background stars, or 
because of small fluctuations in the true background in- 
tensity that are not captured by the large scale averaging 
used in the model. Examining the distribution of intensi- 
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Fig. 2. — Top left: Spitzer IRAC 8/im image of Cloud C (G028. 37+00. 07) (logarithmic intensity scale in units of MJy sr _1 ), with the 
dashed ellipse defined by Simon et al. (2006a) (based on MSX images), the crosses showing mm emission cores (Rathborne et al. 2006), 
and the horizontal line showing a scale of 3'and the intensity scale in MJy sr _1 . Middle left: Large-scale Median Filter (LMF) estimate 
°f A/.o.obs (linear intensity scale in units of MJysr -1 ). Bottom left: LMF estimate of S (intensity scale in units of gcm~ 2 ). Top right: 
A/,o,obs,LMF/-fi/.o,obs,SMF- Middle right: Small-scale Median Filter (SMF) estimate of i"i,,o,obs- Bottom right: SMF estimate of S. 



50 100 150 200 

I M . ob8 (MJy sr-0 

Fig. 3. — Distribution of pixel intensities for Cloud C. The solid 
line shows the distribution inside the LMF (13' by 13') centered on 
the IRDC. The dotted line shows the distribution inside the IRDC 
ellipse and the dashed line shows it outside the ellipse (for a region 
24' square). 





0.8 0.9 1 1.1 1.2 

L,o,lmf/^,o,smf 

Fig. 4. — Distribution of the ratio (/„,o,LMr/^,o,SMF) of back- 
ground intensities derived from of LMF and SMF methods. The 
dotted and dashed lines show the distribution inside and outside 
(for a region 24' square) the cloud cllispsc, respectively. 



ties in the LMF, we find the FWHM of the distribution. 
We set those pixels with l v .\ > 7„,o + 0.5FWHM, equal 
to /„,o so that the derived S at this location is zero. 
These regions are typically stars or bright emission re- 
gions. The remaining pixels with I u _\ > /„,o are allowed 
to yield a negative value of E, which helps prevent small 
scale fluctuations (which can include instrument noise) 
from biasing the total mass in a region to positive values. 

3.2.2. Small-Scale Median Filter (SMF) 

In order to resolve smaller-scale variations in the IR 
background, but without having the median filter es- 
timate be affected by the IRDC, we develop a second 
method, which uses a small-scale median filter in regions 
outside a defined cloud boundary and an interpolation 
scheme inside this boundary. For cloud boundaries, we 
utilize the ellipses defined by Simon et al. (2006a) from 
their study of MSX images of the clouds. We set the filter 
size to be one third of the major axis of this ellipse. The 
method of estimating the median pixel intensity used for 
the LMF (§3.2.1) is repeated for the region outside the 
cloud ellipse. 

For each image pixel inside the IRDC ellipse we esti- 
mate obs by interpolating from values of /„.o,obs out- 
side of the ellipse that are within an angular distance 
equal to the semi-major axis of the ellipse, weighting by 
the inverse square of the angular separation, so that the 
innermost annuli dominate the average. This Small-Scale 
Median Filter (SMF) method is illustrated for cloud C in 
the right-hand panels of Figure 2. Here, and also in Fig- 
ure 4, we compare the LMF and SMF estimates of I v $, 
finding typical variations at the level of < 10%. A 10% 
uncertainty in I v $ corresponds to a mass surface density 



of 0.013 gem -2 . We conclude that systematic uncertain- 
ties in background estimation set a minimum threshold 
E, below which the mid-IR extinction mapping method 
becomes unreliable. 

Here, we also note that the GLIMPSE images are occa- 
sionally prone to certain artifacts, e.g. changes in diffuse 
intensity along diagonal bands and stripes (see cloud B 
and D images below), which will introduce additional un- 
certainties in these regions. The SMF method improves 
upon the LMF method in reducing the effect of these 
artifacts. 

3.2.3. Orthogonal Strips Across Filamentary IRDCs 

To further test the accuracy of the LMF and SMF 
methods, we consider the structure of the IR background 
on very small scales across a filamentary IRDC (cloud 
H). We choose two strips, approximately orthogonal to 
regions of the IRDC that are very filamentary (i.e. thin) 
and free of stars (see Figure 5). The median image inten- 
sities are evaluated along these strips (Figure 6). A lin- 
ear fit for the IR background is fit to regions of the strip 
judged to be free of cloud material. This is a somewhat 
subjective process, since it is hard to distinguish between 
small-scale background variations and additional absorb- 
ing components. Nevertheless, we expect this estimate of 
the background intensity just behind the filament to be 
more accurate than the LMF and SMF methods, because 
the interpolation behind the filament is over a very small 
angular scale (about a few tens of arcseconds) and the 
background appears to be quite smooth in this general 
region. 

The ratio of the LMF and SMF estimates of the back- 
ground across the strips compared to that directly es- 
timated from the strips is shown in Figure 6. We also 
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Fig. 5.— GLIMPSE 8/im image of Cloud H with 3' scale bar 
showing the strips considered in §3.2.3 (logarithmic intensity scale 
in units of MJy sr _1 ). 



show the derived column densities. We see that the vari- 
ations in the background estimates are relatively small, 
at approximately the 10% level, similar to the variations 
between the LMF and SMF methods seen in cloud C, 
and again corresponding to mass surface density uncer- 
tainties of ~ 0.013 g cm~ 2 . 

4. RESULTS 

4.1. Mass Surface Density Maps and Distributions 

The mass surface density maps derived using the SMF 
method for the remaining 9 clouds are presented in Fig- 
ures 7 — 9 with a uniform scale range in E. These IRDCs 
exhibit a variety of morphologies, ranging from very fila- 
mentary (clouds F and H) to those with more apparently 
spherical distributions (clouds C and E). The derived 
mass surface densities range up to ~ 0.35 g cm -2 , which 
is likely to be dependent on the angular resolution of the 
images. As is apparent from the 8 /im images, some of the 
clouds have a few apparently embedded sources, which 
cause localized regions in the clouds to be IR-bright, and 
thus not amenable to our extinction mapping technique. 

In figure 10 we show the mass- weighted probability dis- 
tribution functions (PDFs) with S (evaluated with the 
SMF method) of the 10 IRDCs using the regions inside 
their Simon et al. (2006a) ellipses. Here we define M' 
as the mass fraction, normalized over the distribution of 
pixels with S > 0. In figure 11 we show the same PDF 
but now as a function of E/E, where E is the mean E 
(area- weighted over pixels which have E > 0). The val- 
ues of E are listed in Table 1 . 

The distributions of E/E show the most similarity for 
E/E < 0.5. The high-side of the distributions do show 
a wide variety of profiles: there is a large range of mass 



fractions at mass surface densities that are, for example, 
3 times greater than the mean. Cloud J has the narrowest 
PDF, while cloud C has the broadest and one of the 
most skewed distributions. We caution that our derived 
IRDC PDFs are likely to be sensitive to the geometry of 
the ellipse that was chosen to represent the lRDC. For 
example this can affect the normalization of E (although 
note we restrict its calculation over only the area of the 
cloud for which we measure E > 0). 

In the lower panel of Figure 11 we make a qualitative 
comparison of these results with the shapes of E/E PDFs 
resulting from different astrophysical sources of turbu- 
lence as modeled in two sets of published numerical simu- 
lations. Offner, Klein, & McKee (2008) presented simula- 
tions of isothermal driven turbulence to maintain a mean 
ID Mach number of about 5 in which self-gravity was 
then turned on at a particular time. One expects clouds 
with higher Mach number turbulence to have broader 
PDFs of E/E, as stronger shocks lead to greater con- 
trasts between regions of compression and rarefaction. 
Clouds that are more strongly self-gravitating also have 
broader PDFs. Self-gravity also tends to skew the PDFs 
towards the high- value side. Figure 11 shows two exam- 
ples from these simulations: one is before self-gravity has 
been initiated, so is representative of the E/E PDF from 
hydrodynamic Mach 4.7 (ID) driven turbulence with a 
Burgers P{k) cx k~ 2 power spectrum; the other, which 
has a similar Mach number of 4.9, is after self-gravity 
has been allowed to develop structures for about 1 mean 
simulation free-fall time while the turbulence is still be- 
ing driven and is somewhat broader as a result (14.2% 
of the gas mass has collapsed into sink particles by this 
time; the PDF is normalized to the total remaining gas 
mass). 

Nakamura & Li (2007) presented simulations of mag- 
netized protostellar outflow driven turbulence. The driv- 
ing scale is relatively small compared to the simulation 
box, and the sources are centrally concentrated in the 
box. Their formation is sporadic as the global clump 
collapses: we show results after 1.5 global gravitational 
collapse times, when about 80 protostars have formed. 
Furthermore there is a large scale direction to the initial 
magnetic field, which has a mean dimensionless flux-to- 
mass ratio of 0.52 (in units of the critical value, 2ttG x ^ 2 
[Nakano & Nakamura 1978]), decreasing to 0.19 in the 
clump center. The mean value corresponds to about 
75 fiG for the fiducial simulation parameters. This in- 
hibits perpendicular motions, as can be seen from the 
narrowness of the E/E PDF derived looking along the 
field direction compared to perpendicular to it (Fig- 
ure 11; these distributions correspond to the models 
shown in Figure 7 of Nakamura & Li [2007]). 

A comparison of the observed and simulated E/E 
PDFs shows that the ensemble of the observed distri- 
butions can be qualitatively accounted for by the range 
of simulations shown, although the observed IRDCs gen- 
erally have narrower distributions. However, we arc not 
able to exclude the possibility that other numerical mod- 
els involving different physics, e.g. a larger degree of 
magnetic support, and/or other parameters could not 
also explain the observed clouds. Since most of the 
IRDCs have relatively narrow PDFs compared to the 
Mach 5 turbulence models, this may indicate that these 
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Fig. 6. — Results of LMF and SMF methods with thin strip analyses for strip a (left column) and strip b (right column). Top panels: 
Median i b s along the strip (solid lines), LMF estimate of -f^o.obs (dotted lines), SMF estimate of o,obs (dashed lines), estimate from 
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LMF (dotted) and SMF (dashed) estimates of /^o.obs compared to estimate from the strip intensity profiles. Bottom panels: Estimates of 
S via LMF (dotted), SMF (dashed), and strip analysis (long dashed) methods. 



clouds have smaller turbulent velocity dispersions and/or 
that dynamically important magnetic fields are present. 
The protostellar outflow driven turbulence models gener- 
ally have broader S/S PDFs than the observed IRDCs, 
and this may indicate these models are more applica- 
ble to later evolutionary stages of star cluster formation. 
However, they do illustrate the effects of dynamically im- 
portant large-scale magnetic fields in creating narrower 
S/S PDFs for viewing directions along the field lines. 

Observations of the velocity structure of these IRDCs 
as traced by molecular line emission can help to further 
test the nature of turbulence present in the early stages 
of star cluster formation. The extinction mapping tech- 
nique we have presented also opens up the possibility 
of studying the E/S PDF as a function of IRDC £, M 



and embedded stellar content. We defer a more quan- 
titative study of these dependencies and a quantitative 
comparison of the S/S PDFs to those formed in numer- 
ical simulations to a future study. 

4.2. Cloud and Core Masses 

Given the (kinematic) distance to each cloud (Simon 
et al. 2006a), we convert our S estimates into mass es- 
timates. For a typical distance of 3 kpc, each 2.0" res- 
olution element at 8 /im (i.e. the mean FWHM of the 
point response function) corresponds to a linear scale of 
0.029 pc, and if this has a typical S of 0.1 g cm~ 2 , 
this corresponds to a mass of 0.41 M Q . We see that the 
mid-IR extinction mapping technique has the potential 
to probe relatively low mass scales. 

In Table 1 we list the total cloud masses derived by 



Fig. 7. — The left column shows Spitzcr IRAC 8fim images of IRDCs, with the dashed ellipse defined by Simon ct al. (2006a) (based on 
MSX images), the crosses showing mm emission cores (Rathborne et al. 2006), with the numbers referring to their designation in Table 3, 
and the horizontal line showing a scale of 3'. The right column shows extinction maps made by the small-scale median filtering method 
(SMF). From top to bottom, the clouds are: Cloud A (G018.82-00.28), Cloud B (G019.27+00.07), Cloud D (G028.53-00.25). 
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Fig. 8.— Same as Figure 7, but now, from top to bottom, the clouds are: Cloud E (G028.67+00.13), Cloud F (G034.43+00.24), Cloud 
G (G034.77-00.55). 
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Fig. 9.— Same as Figure 7, but now, from top to bottom, the clouds are: Cloud H (G035.39-00.33), Cloud I (G038.95-00.47), Cloud J 
(G053. 11+00.05). 



the LMF and SMF extinction mapping methods. Note, 
these are the masses inside the elliptical regions defined 
by Simon et al. (2006a) from MSX images, and thus are 
not necessarily particularly close representations of the 
morphologies revealed by Spitzer IRAC. We find a mean 



fractional difference of 35% between these two mass es- 
timates. We regard the SMF mass estimate as being 
more accurate. Uncertainties in both methods will grow 
for larger clouds with lower mean mass surface densities 
(i.e. smaller contrasts against the background). Note 
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E (thin dot-dashed black), F (thick solid gray [pink in electronic 
version]), G (thick dotted gray), H (thick dashed gray), I (thick 
long-dashed gray), J (thick dot-dashed gray). 



also this extinction mapping technique is not sensitive 
to a uniform screen of matter that covers both the re- 
gion of the IRDC and the region where the background 
is estimated. A comparison of extinction-derived cloud 
masses and mass surface densities with the properties in- 
ferred by molecular line emission, such as 13 CO, will be 
presented in a separate study (A. K. Hernandez & J. C. 
Tan, in prep.). 

The extinction mass estimates are expected to become 
more accurate at high values of S. Each of the IRDCs 
in our sample contains dense cores studied by their mm 
dust continuum emission by Rathborne et al. (2006). We 
first identify those cores which are amenable to extinc- 
tion mapping (i.e. do not contain bright 8 /im emission). 
Cores identified by Rathborne et al. (2006) as being mid- 
IR bright (labelled "(e)" in Table 3) were excluded. We 
also excluded cores overlapping with fainter IR sources, 
if the area affected by the sources was greater than about 
10% of the core. Finally, we excluded cores with very low 
surface densities, £ < 0.02 gem -2 , since the extinction 
mapping method becomes very unreliable at these levels 
due to uncertainties in the background estimation. 

For our selected cores, we summed the mass inside a 
circular area centered on each core position with radius 
equal to half the angular FWHM diameter of Gaussian 
fits that Rathborne et al. (2006) fitted to their mm con- 
tinuum images. Note the quoted angular FWHM diame- 
ters are typically about a factor of 1.09 times larger than 
the deconvolved FWHM diameters, so our areal integra- 
tion actually extends out to about 0.545 true FWHM 
diameters from the core center. The total core masses 
quoted by Rathborne et al., which are evaulated by in- 
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FlG._ll. — Top: Mass-weighted probability distribution function 
of S/S for clouds A (thin solid black line), B (thin dotted black), C 
(thin dashed black), D (thin long-dashed black), E (thin dot-dashed 
black), F (thick solid gray [pink in electronic version]), G (thick 
dotted gray), H (thick dashed gray), I (thick long-dashed gray), 
J (thick dot-dashed gray). Note, cloud J has the most sharply- 
peaked, narrowest distribution, while cloud C has the broadest. 
Bottom: Comparison to simulations of astrophysical turbulence. 
Dotted line shows the result of Offner et al. (2008) for driven 
(Mach 4.7) turbulence with no self-gravity. Solid line shows the 
same simulation after self-gravity has been allowed to operate for 
1 mean free fall time. The long-dashed line shows the result of 
Nakamura & Li (2007) for protostellar driven turbulence after 1.5 
global gravitational collapse times as an average of two orthogonal 
viewing angles that are both perpendicular to the mean large-scale 
B-ficld in the simulation. The dashed line shows the view along 
the mean B-field. 



tegrating over the whole 2D Gaussian distribution, will 
thus be a factor of 1.78 larger than that contained in 
this angular area, so we reduce their masses by a fac- 
tor of 0.562. Their mass estimates involved assuming 
T = 15 K, a dust opacity per mass of dust of Ki.2mm = 
1.0 cm 2 g -1 (Ossenkopf & Henning 1994), a gas-to-dust 
mass ratio of 100, and graybody emission with emissivity 
index (3 = 2. In the context of the Ossenkopf & Henning 
(1994) dust model with thin ice mantles, coagulated for 
10 5 yr at 7Jh = 10 6 cm -3 , which is the closest to our 
adopted 8 /Ltm opacity, the interpolated 1.2 mm opacity 
is 1.056cm 2 g _1 , i.e. 5.6% higher than assumed by Rath- 
borne et al. (2006), so we further reduce their masses by 
a factor 0.947. Our calculation of S from dust extinc- 
tion assumed a gas-to-dust ratio of 156, which was that 
adopted for the Ossenkopf & Henning (1994) coagulation 
model. If this same ratio is applied to the mm emission 
masses, then the derived masses would be raised by 1.56 
from the Rathborne et al. estimates. These three effects 
imply the Rathborne et al. quoted core masses should be 
reduced by a factor of 0.830 for a fair comparison with 
our extinction masses. Note, we use the same distance 
adopted by Rathborne et al. to the cores. In Table 3 and 
Figure 12 we compare the mm dust emission and mid-IR 
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dust extinction mass estimates. 

We find generally very good agreement between these 
different methods. The dispersion in the ratios of the 
LMF and SMF estimates is about 15%. Comparing 
to the mm emission masses, the logarithmic average of 
M mm /MsMF is 1.08 with a dispersion of a factor of 1.8. 
This small systematic offset may be somewhat fortuitous 
given the systematic errors inherent to both methods. 
The possibility of correlations of M mm /M SMF with core 
properties is discussed below. 

4.3. Correlations of Core M mm /MsMF with Density as 
Evidence for Grain Growth 

We also use our SMF extinction mass to calculate 
a volume- aver aged H number density in each core, 
^h,smf, assuming spherical geometry. We find values 
in the range nn.SMF — 10 4 — 10 5 cm~ 3 . We exam- 
ine the ratio of M mm /MsMF as a function of core den- 
sity in Figure 13. The best fit power law relation is 
(M mm /MsMF) oc Wjj° s ^ F . Considering the Spearman 
rank-order correlation, the probability of a chance corre- 
lation is 0.0036. Since Msmf and tjh,smf are correlated 
via the observable Xsmf, we also examine the correlation 
of (M mm /M SMF ) with M S mf, Esmf (Figure 12), core ra- 
dius, r, and core distance (Figure 13). We find proba- 
bilities of chance (anti)correlation of 0.28, 0.018, 0.99, 
and 0.43 respectively, i.e. there is no significant cor- 
relation of (M mm /MsMp) with Msmf, r, and distance, 
and there is a marginally significant anticorrelation with 
Ssmf- The most significant trend is the anticorrelation 
of (M mm /MsMF) with density. 

One possible explanation for this observed trend is 
a systematic change in dust opacities at 8fim and 
1.2 mm caused by grain growth, since M mm /MsMF oc 
%m/Ki.2mm- For example, the thin ice mantle model of 
Ossenkopf & Henning (1994) predicts KiRAC8/im/Ki.2mm 
decreases by a factor of 0.68 after coagulation for approx- 
imately 10 6 yr at a density of 71h = 10 5 cm~ 3 (or 10 7 yr 
at a density of nn = 10 4 cm~ 3 ). The ratio decreases only 
slightly in the limit of further coagulation. Compari- 
son of the uncoagulated bare grain and thin ice man- 
tle models of Ossenkopf & Henning (1994) shows that 
K8jum/Ki.2mm decreases by about a factor of 0.75 by the 
formation of ice mantles. Thus a total factor of about 
0.5 decrease in M mm /MsMF, i.e. about that observed in 
our cores as the mean density changes from ~ 10 4 cm~ 3 
to ~ 10 cm , could be explained by a combination of 
ice mantle formation and grain growth. Note the typi- 
cal density in the cores could be higher than the mean 
volume averaged value if they are clumpy. 

Systematic temperature changes may also lead to 
changes in M mm /MsMF, via the estimate of M mm . One 
does expect colder temperatures at higher densities, 
which would lead to M mm being relatively underesti- 
mated at high densities under the constant temperature 
assumption of Rathborne et al. (2006). For example, 
adopting a temperature of 10 K rather than 15 K would 
raise M mm by a factor of 1.9, which is enough to account 
for the observed trend of M mm /MsMF with density. 

We conclude that this IRDC core population, as the 
mean density increases from ~ 10 4 cm~ 3 to ~ 10 5 cm~ 3 , 
shows tentative evidence for either opacity changes that 
are consistent with models of ice mantle formation and 



grain coagulation, or a temperature decrease of ~ 5 K. 
Systematic follow up of these cores with NH3 observa- 
tions to measure temperature, as has been done in some 
IRDCs (e.g. Sridharan et al. 2005; Pillai et al. 2006, 
Wang et al. 2008), can help to distinguish these possi- 
bilities. 

Evidence for grain growth has previously been reported 
in the low- mass core B68 by Bergin et al. (2006), based 
on low gas-to-dust temperature coupling rates. Bianchi 
et al. (2003) considered ratios of 850 /im and 1.2 mm 
opacities to visual opacities in B68, finding no conclusive 
evidence for grain growth. Keto & Casclli (2008) require 
increased dust opacities due to large, fluffy grains to ex- 
plain the < 7 K temperatures observed in the centers of 
some low-mass starless cores. Flower, Pineau des Forets, 
& Walmsley (2005, 2006) considered the chemical impli- 
cations of grain growth, finding it causes an increase in 
the electron fraction and the ratio of H + /H^, and can 
have moderate effects on the observed differential freeze- 
out of nitrogen and carbon bearing species. Vastel et 
al. (2006) required grains larger than 0.3^m to repro- 
duce the ortho-H 2 D+ observations in the prestcllar core 
L1544. A systematic molecular line study of our studied 
IRDC cores to test for the above chemical effects and any 
correlation with M mm /MsMF would be useful. 

5. CONCLUSIONS 

We have presented a new method of diffuse mid-IR 
(8/im) extinction mapping of infrared dark clouds to 
probe the initial conditions for massive star and star 
cluster formation. The technique can naturally probe to 
higher mass surface densities than NIR extinction map- 
ping using background stars and provides an estimate 
of the mass surface density at scales that are only lim- 
ited by the resolution of the mid-IR image (e.g. a couple 
of arcseconds for Spitzer IRAC). However, the method 
involves a spatial interpolation and smoothing of the dif- 
fuse background, which inevitably introduces additional 
uncertainties. In particular it is a relative mass sur- 
face density of the cloud compared to the region used 
to ascertain the background that is derived. Also, small 
scale background fluctuations are not captured and lo- 
cal regions containing mid-IR bright sources cannot be 
treated. One of the largest uncertainties for more distant 
clouds is the correction for foreground mid-IR emission 
from hot dust, and so we have generally restricted the 
analysis to relatively nearby clouds. 

We have examined three different methods of estimat- 
ing the diffuse background behind IRDCs. The Large- 
scale Median Filter (LMF) method is the simplest by not 
requiring information about the cloud size or location, 
but then must smooth over relatively large scales. The 
Small-scale Median Filter (SMF) achieves higher spatial 
resolution of the background given defined cloud bound- 
aries, and we generally regard this method as being more 
accurate where such information is available. Finally, for 
very thin, filamentary clouds, more detailed background 
fitting can be attempted by considering perpendicular 
strips across the filament. Comparing these methods on 
a sample of 10 IRDCs we find systematic differences due 
to background fitting at about the 10% level. This means 
that the mid-IR extinction mapping technique becomes 
unreliable for E < 0.013 g cm~ 2 (given our adopted fidu- 
cial 8 fim [Spitzer IRAC Band 4 and diffuse Galactic 
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Fig. 12. — Comparison of core mass estimates. Top Left: Ratio of cores masses determined via extinction mapping with Large and Small 
Median Filter methods (Mlmf/Msmf) versus core mass. Top right: Same data plotted as a function of core S (estimated via the SMF 
method). Middle panels: Ratio of mm emission mass (Rathborne et al. 2006) to Mlmf- Bottom panels: Ratio of mm emission mass to 
M SMF . 



background spectrum weighted] opacity of 7.5 cm 2 g _1 ). 

We have then used this method to measure the mass 
surface densities of the 10 IRDCs. The probability dis- 
tribution functions of S for these clouds show a range 
of shapes, that is qualitatively consistent with that ex- 
pected in numerical simulations of turbulence, though 
the relatively narrow distributions of the observed clouds 
indicate that the Mach numbers are probably < 5 and/or 
that magnetic fields are dynamically important. A more 
quantitative comparison of the observed clouds with nu- 
merical models, and a search for evolutionary trends as- 
sociated with the development of self-gravitating struc- 
tures is planned for a future study. 

Given IRDC kinematic distances, we then estimate 
masses of the clouds and, more accurately, of their dense 
cores. Comparing to mass estimates from mm dust con- 
tinuum emission (Rathborne et al. 2006) we find very 



good agreement: a < 10% systematic offset for the pop- 
ulation of 43 studied cores, and a dispersion of less than 
a factor of two in the distribution. 

Finally, we examine the ratio of masses estimated by 
mm dust emission to those found from mid-IR dust ex- 
tinction as a function of core density, tentatively finding 
a trend which can be explained as being due to opacity 
changes due to ice mantle formation and grain growth 
or by a temperature decrease of about 5 K in the dens- 
est cores. Future studies of the mid-IR extinction law in 
these cores and molecular line studies of their chemistry 
can help to test the reality of this result and distinguish 
between these explanations. 
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Fig. 13. — Left panel: Ratio of core masses estimated from mm dust emission and 8 ^tm extinction with the SMF method (M mm /MsMF) 
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of this anticorrelation is discussed in the text. Middle panel: Ratio of (M mm /MsMF) versus core radius. No significant trend is observed 
(see text). Right panel: Ratio of (M mm /MsMF) versus core distance. No significant trend is observed (see text). 
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TABLE 3 
IRDC Cores 



Core Namc a 


Galactic Coordinates 


Diamctcr b 


M mm c 


^LMF 


SsMF 




M SMF 


"H,SMF 




/ 


b 


(PC) 




(g cm -2 


) (gem- 2 ) 


(Mq) 




(10 s cm- 3 ) 


Al (MM4) 


18.790 


—0.286 


0.721 


219 


0.0829 


0.0813 


163 


157 


0.232 


A2 (MM6) 


18.799 


—0.294 


0.814 


171 


0.0805 


0.0710 


211 


188 


0.193 


A3 (MM5) 


18.806 


—0.303 


0.698 


141 


0.0805 


0.0749 


156 


142 


0.230 


A4 (MM3)(e) 


18.701 


—0.227 


0.582 


194 


- 


- 


- 


- 


- 


A5 (MMlj(c) 


18.735 


—0.226 


0.535 


736 


- 


- 


- 


- 


- 


A6 (MM2)(c) 


18.833 


—0.299 


0.465 


201 


- 


- 


- 


- 


- 


Bl (MM2) 


19.288 


0.0824 


0.431 


94.5 


0.0656 


0.0631 


57.3 


51.3 


0.354 


B2 (MM1) 


19.311 


0.0675 


0.337 


93.7 


0.0774 


0.0687 


32.2 


25.7 


0.370 


CI (MM9) 


28 324 


n nfi77 


0.824 


330 


0.165 


0.194 


401 


539 


0.532 


C2 (MM4) 




u.uoy ( 


0.509 


274 


0.234 


0.304 


239 


303 


1.26 


C3 (MM11) 




n 1 nns 

U. 1UUO 


0.921 


371 


0.199 


0.235 


243 


279 


0.197 


C4 (MM6) 


ZO.OJO 


U.U I 


0.485 


192 


0.162 


0.186 


204 


248 


1.20 


C5 (MM14) 


OQ qcifi 


u.uouu 


0.291 


30.0 


0.188 


0.238 


66.1 


71.1 


1.59 


C6 (MM10) 


28.362 


0532 


0.824 


296 


0.0757 


0.0857 


465 


584 


0.576 


C7 (MM16) 


OQ ocj 


U. -LZUO 


1.07 


371 


0.203 


0.256 


352 


406 


0.183 


C8 (MM17) 


OQ Q Q Q 
ZO.OOO 


U.UOOl 


0.727 


164 


0.0789 


0.0892 


329 


387 


0.556 


C9 (MMl)(e) 


OQ QQQ 


n nsi 9 


0.630 


952 


0.155 


0.183 


248 


291 


0.643 


CIO (MM8)(e) 


OQ OAA 


n ni 97 


0.751 


343 


- 


- 


- 


- 


- 


Cll (MM7)(e) 


OQ OQO 


u.uuuo 


0.582 


253 


- 


- 


- 


- 


- 


C12 (MM3)(e) 


OQ Q99 


— U.UlUi 


0.558 


400 


- 


- 


- 


- 


- 


C13 (MM5)(e) 


OQ Q9^ 


U. -LUi-O 


0.412 


147 


- 


- 


- 


- 


- 


C14 (MM12)(e) 


OQ QOQ 


— U.UoOo 


1.02 


395 


- 


- 


- 


- 


- 


C15 (MM2)(e) 


OQ QQ7 
ZO.DO 1 


U. J. 1 l u 


0.533 


449 


- 


- 


- 


- 


- 


C16 (MM15)(e) 


28.339 


0.1424 


0.630 


111 


- 


- 


- 


- 


- 


C17 (MM18)(e) 


28.417 




0.824 


135 


- 


- 


- 


- 


- 


CI8 (MM13)(c) 


28.419 


0.1391 


1.02 


385 


■ 


- 




- 


■ 


Dl (MM5) 


OQ ^9fi 




0.636 


194 


0.0873 


0.0995 




126 


0.271 


D2 (MM7) 


9°. ^8 


—0.2757 


1.24 


421 


0.0799 


0.0942 


. 


441 


0.127 


D3 (MM3) 


OS C^/1Q 


9^fiQ 

— u.zouy 


1.55 


1690 


0.0722 


0.0844 


672 


729 


0.108 


D4 (MM8) 


28.543 


—0.2651 


0.746 


154 


0.0565 


0.0669 


119 


127 


0.169 


D5 (MM2) 


28.559 


-0.2274 


1.55 


1760 


0.0817 


0.0924 


753 


813 


0.121 


D6 (MM4) 




—0.2412 


1.22 


663 


0.0707 


0.0739 


324 


380 


0.116 


T~""\ r~7 /T\ TT\ T1 \ 

D7 (MM1) 


28 565 




0.912 


1000 


0.0603 


0.0655 


280 


318 


0.232 


D8 (MM10) 


28.579 


—0 2303 


1.13 


360 


0.0708 


0.0836 


435 


447 


0.171 


D9 (MM9) 


28 589 


—0 2285 


0.718 


138 


0.0829 


0.0959 


140 


160 


0.239 


D10 (MM6) 


28.557 


—0 9^89 


0.497 


98.6 


- 


- 


- 


- 


- 


El (MM7) 


OQ fi/l/L 


fl 1 ^7^ 


0.742 


72.2 


0.109 


0.119 


206 


230 


0.311 


E2 (MM5) 


OQ fi^n 


fl 1 9fifl 


0.742 


85.4 


0.0632 


0.0774 


193 


221 


0.298 


E3 (MM2) 


OQ fifil 
Zo.OOl 


fl 1 A^fi 


0.940 


327 


0.0872 


0.100 


359 


398 


0.264 


E4 (MM4) 


28.717 


0.1456 


0.717 


96.2 


0.0953 


0.106 


140 


169 


0.254 


E5 (MM3) 




0.1456 


0.544 


101 


- 


- 


- 


- 


- 


E6 (MM6) 


28.647 


0.1143 


0.766 


93.7 












E7 (MM1) 


9k fi&S 
zo.uoo 


0.1782 


0.495 


119 


- 


- 




- 


- 


Fl (MM8) 


34.422 


0.24792 


0.556 


89.6 


0.0556 


0.0497 


- 


67.8 


0.217 


F2 (MM7) 


34.438 




0.502 


72.2 


0.0573 


0.0517 


55.8 


49.8 


0.217 


b 3 [MMbJ 


34.448 


2 r S0Q1 


0.664 


104 


0.0625 


0.0591 


97.0 


88.4 


0.166 


F4 (MM9) 


34.454 


fl 25495 


0.699 


130 


0.0583 


0.0513 


116 


89.7 


0.145 


F5 (MM5) 


34.386 


21 378 


0.915 


550 


- 


- 


- 


- 


- 


F6 (MM4) 


34.394 


n 221 56 


0.431 


211 


- 


- 


- 


- 


- 


F7 (MM2) 


34.401 


n 22709 

KJ . ZiZi 1 KJis 


0.466 


1060 


- 


- 


- 


- 


- 


F8 (MM1) 


34.411 


23489 


0.287 


985 


- 


- 


- 


- 


- 


F9 (MM3) 


34.459 


24867 


0.431 


250 


- 




- 


- 


- 


Gl (MM3) 


34.734 


—0.5670 


0.225 


11.6 


0.0606 




10.4 


9.77 


0.474 


G2 (MM2) 


34.783 


—0 ^683 

•J. UUOJ 


0.689 


156 


0.0515 


■ 


112 


105 


0.177 


C3 (MM4) 


34.784 


—0 ^608 


0.506 


58.0 


0.0667 


0.0682 


64.1 


64.6 


0.275 


G4 (MM1) 


34.712 


—0.5946 


0.323 


138 


- 


- 


- 


- 


- 


HI (MM9) 


35.478 


-0.3096 


0.380 


34.8 


0.0842 


0.0869 


78.7 


69.9 


0.705 


H2 (MM4) 


35.483 


-0.2858 


0.478 


89.6 


0.113 


0.114 


69.9 


71.9 


0.364 


rlo ^IVilVlO ) 


35.483 


-0.2954 


U.ozU 


y / .y 


U.UDUz 


U.UDZ1 


i in 

iiy 


1 OO 


U.4 it 


H4 (MM8) 


35.491 


-0.2830 


0.408 


48.9 


0.0830 


0.0854 


36.0 


37.4 


0.304 


H5 (MM6) 


35.497 


-0.2863 


0.520 


58.8 


0.0602 


0.0629 


57.8 


65.3 


0.257 


H6 (MM7) 


35.522 


-0.2724 


0.548 


79.6 


0.133 


0.138 


89.8 


92.0 


0.308 


H7 (MM2)(c) 


35.417 


-0.2847 


0.239 


37.3 












H8 (MM3)(c) 


35.452 


-0.2951 


0.337 


65.5 












H9 (MMl)(c) 


35.456 


-0.1801 


0.295 


63.0 












11 (MM1) 


38.957 


-0.4659 


0.380 


97.0 


0.111 


0.124 


57.3 


63.6 


0.640 


12 (MM3) 


38.971 


-0.4588 


0.157 


9.12 


0.135 


0.142 


12.8 


11.6 


1.65 


13 (MM2)(e) 


38.937 


-0.4578 


0.327 


60.5 












14 (MM4)(c) 


38.949 


-0.4385 


0.432 


39.8 
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TABLE 3 — Continued 



Core Namc a 


Galactic Coordinates 


Diamctcr b 


M mm c 


Slmf Ssmf 


M-LMF MsmF 


«H,SMF 




/ 


b 


(pc) 


(M e ) 


(gem" 2 ) (gem" 2 ) 


(M e ) (M e ) 


(10 5 cm- 3 ) 


Jl (MM4) 


53.128 


0.0503 


0.384 


37.3 


0.0324 0.0219 


17.1 12.2 


0.118 


J2 (MM3)(e) 


53.092 


0.1201 


0.166 


9.94 








J3 (MM5)(e) 


53.136 


0.0283 


0.209 


10.8 








J4 (MMl)(c) 


53.141 


0.0717 


0.183 


103 








J5 (MM2)(e) 


53.157 


0.0672 


0.288 


36.5 









a Name in parentheses from Rathborne ct al. (2006). "(c)" indicates their designation of the core as containing 8 fim emission. We have 
identified additional cores that contain 8 fim sources, for which we also do not calculate extinction mass surface densities and masses. 
b Derived from the angular size of the FWHM of the Gaussian profile fitted to the core by Rathborne et al. (2006). 
c Millimeter emission mass, from Rathborne ct al. (2006). 



